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Abstract. We present a new evaluation of an Sr optical lattice clock using spin 
polarized atoms. The frequency of the ^So ^Po clock transition is found to be 
429 228 004 229 873.6 Hz with a fractional accuracy of 2.6 x 10"^^ a value that 
is comparable to the frequency difference between the various primary standards 
throughout the world. This measurement is in excellent agreement with a previous 
one of similar accuracy [1] . 



1 Introduction 

The possibility to build high accuracy optical clocks with ^''Sr atoms confined in an optical 
lattice is now well established. Since this idea was published [2] , experiments rapidly proved 
the possibility to obtain narrower and narrower resonances with atoms in the Lamb-Dicke 
regime |3l4l5j . The narrowest observed resonances now have a width in the Hz range [6j and 
the corresponding potential fractional frequency instabilities are better than 10""'^^ over 1 s of 
averaging time. On the other hand, systematic effects were also shown to be highly control- 
lable. It was theoretically demonstrated that the residual effects of the atomic motion could 
be reduced down to the 10"^* level for a lattice depth as small as lOE'r with Er the recoil 
energy associated with the absorption or emission of a lattice photon. Higher order frequency 
shifts due to the trapping light were then also shown to be controllable at that level [5pl . Alto- 
gether, the accuracy of the frequency measurement of the ^Sq — ^Pq clock transition of Sr has 
steadily improved by four orders of magnitude since its first direct measurement in 2003 [5] . 
Three independent measurements performed in Tokyo university [9] , JILA [4] and SYRTE |10j 
were reported with a fractional uncertainty of 10^^** giving excellent agreement. Recently, the 
JILA group improved their uncertainty down to 2.5 x 10~^^ [T], a value that is comparable to 
the frequency difference between the various primary standards throughout the world |llj . We 
report here a new and independent measurement of this clock transition with an accuracy of 
2.6 X 10~^^. The major modification as compared to our previous evaluation is an improved 
control of the Zeeman effect. By applying a bias field of typically 0.1 mT and pumping atoms 
into extreme Zeeman states (we alternate measurements with mp = +9/2 and niF = —9/2) 
we cancel the first order Zeeman effect while getting a real time measurement of the actual 
magnetic field seen by the atoms 'W . The measured frequency of the ^Sq — > '^Pq clock transition 
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of s^Sr is 429 228 004 229 873.6 (1.1) Hz. This value differs from the one of Ref. T by 0.4Hz 
only. 



2 Experimental setup 
2.1 Atom manipulation 

The apparatus is derived from the one described in Ref. [TD] . The clock is operated sequentially, 
with a typical cycle duration of 400 ms. We use a dipole trap formed by a vertical standing wave 
at 813.428 nm inside an enhancement Fabry-Perot cavity. The depth of the wells of the resulting 
lattice is typically 100 fiK. A beam of ®^Sr atoms from an oven at a temperature of about 450 °C 
is sent through a Zeeman slower and then loaded into a magneto-optical trap (MOT) based on 
the — > ^Pi transition at 461 nm. The MOT temperature is about ImK. The dipole trap 
laser is aligned in order to cross the center of the MOT. Two additional laser beams tuned to the 
^Sq '^Pi and ^Pi '^Si transitions, at 689 nm and 688 nm respectively, are superimposed on 
the trapping laser. The atoms that cross these beams are therefore drained into the metastable 
states '^Pq and '^P2 at the center of the trap, and those with a small enough kinetic energy 
remain confined in the potential wells forming the lattice. The MOT and drain lasers are then 
switched off and atoms are optically pumped back to the ground state, where they are further 
cooled using the narrow ^Sq ^ ^Pi transition. About 95% of the atoms are cooled down 
to the ground state of the trap (see Fig.[T]), corresponding to a temperature of 3 /zK. They are 
optically pumped into either the (^5*0, mp = 9/2) or {^Sq, mp — —9/2) Zeeman sub-state. This 
is achieved by means of a bias magnetic field of about = 10~^ T and a circularly polarized laser 
(cr^ or a~ depending on the desired mp state) tuned to the ^So{F = 9/2) ^ ^Pi{F = 9/2) 
transition. This transition is power-broadened to a few hundreds of kHz. The magnetic field 
can then be switched to a different value (up to a fraction of a mT) for the clock transition 
interrogation. We use a 7r-polarized laser at 698 nm to probe the "^Pq transition with 

adjustable frequency to match the desired (m^ = ±9/2 itif = ±9/2) transition. Finally the 
populations of the two states ^Sq and ^Pq are measured by laser induced fluorescence using two 
blue pulses at 461 nm separated by a repumping pulse. 



2.2 Measurement scheme 

The spectroscopy of the clock transition is performed with an extended-cavity diode laser at 
698 nm which is pre-stabilized with an interference filter [T^] . The laser is stabilized to an ultra- 
stable cavity of finesse F = 25000, and its frequency is constantly measured by means of a 
femtosecond fiber laser |13|14j referenced to the three atomic fountain clocks FOl, F02 and 
FOM. The femtosecond fiber laser is described in paragraph l2.3l 

The fountain ensemble used in this measurement is described extensively in [15|16j . The 
three atomic fountains FOl, F02 and FOM are used as primary frequency standards measuring 
the frequency of the same ultra-low noise reference derived from a cryogenic sapphire oscillator. 
Practically, the reference signal at 11.98 GHz is divided to generate a 100 MHz reference which 
is disseminated to FOl and FOM. Being located in the neighboring lab, F02 benefits from using 
the 11.98 GHz directly. Another 1 GHz reference is also generated from the 11.98 GHz signal 
and sent through a fiber link as a reference for the fiber femtosecond optical frequency comb. 
The 100 MHz signal is also compared to the 100 MHz output of a H-maser. A slow phase-locked 
loop (time constant of 1000 s) is implemented to ensure coherence between the reference signals 
and the H-maser to avoid long term frequency drift of the reference signals. 

During the 15 days of measurement reported here, the three fountains are operated con- 
tinuously as primary frequency standards measuring the same reference oscillator. The overall 
frequency instability for this measurement is 3.5 x 10"^'* at Is for F02, 4.2 x 10"^'' at Is for 
FOl and 7.2 x 10"" at Is for FOM. The accuracy of these clocks are 4 x lO'^^ for FOl and 
F02 and 1.2 x 10^^^ for FOM. The fractional frequency differences between the fountain clocks 
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Fig. 1. Spectrum at high power of the carrier and the first two longitudinal sidebands of the trapped 
atoms. The ratio between both sidebands is related to the population of the ground state of the trap. 
95 % of the atoms are in the lowest vibrational state of the lattice wells. 



are all consistent with zero within the combined 1-sigma error bar, which implies consistency 
to better than 10"^^. 



The link between the cavity stabilized laser and the frequency comb is a fiber link of 50 m 
length with a phase noise cancellation system. The probe laser beam is sent to the atoms after 
passing through an acousto-optic modulator (AOM) driven at a computer controlled frequency. 
Each transition is measured for 32 cycles before switching to the transition involving opposite 
mp states. Two digital servo-loops to both atomic resonances therefore run in parallel with 
interlaced measurements. For each servo-loop we alternately probe both sides of the resonance 
peak. The difference between two successive transition probability measurements constitutes 
the error signal used to servo-control the AOM frequency to the atomic transition. In addition, 
we interlace sets of 64 cycles involving two different trapping depths. The whole sequence is 
repeated for up to one hour. 

This operating mode allows the independent evaluation of three clock parameters. The 
difference between the frequency measurements made for each Zeeman sub-state can be used to 
accurately determine the magnetic field. As we switch to the other resonance every 32 cycles, 
this gives a real-time calibration of the magnetic-field averaged over 64 cycles. The global 
average of the measurement is the value of the clock frequency and is independent on the first 
order Zeeman effect as the two probed transitions are symmetrically shifted. Finally, the two 
frequencies corresponding to two different dipole trap depths are used for a real-time monitoring 
of the possible residual light shift of the clock transition by the optical lattice. 

The frequency stability of the Sr lattice clock-F02 fountain comparison is shown in Fig. [21 
The Allan deviation is 6 x 10~^^t~^/^ so that the statistical uncertainty after one hour of 
averaging time is 10^^^, corresponding approximately to 0.5 Hz. 
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Fig. 2. Allan standard deviation of the frequency measurements for a magnetic field B = 0.87 G and 
a time of interrogation of 20 ms. The line is a fit to the data using a r"^''^ law. The corresponding 
stability at Is is 6 x 10~^*. 



2.3 The frequency comb 

For the absolute frequency measurement of the Sr transition we have used a fibre-based optical 
frequency comb which is based on a FC1500 optical frequency synthesizer supplied by Menlo 
Systems. The laser source for the FC1500 comb is a passively mode-locked femtosecond fibre 
laser which operates at a centre wavelength of approximately 1550 nm and has a repetition rate 
of 100 MHz. The repetition rate can be tuned over approximately 400kHz, by means of an end 
mirror mounted on a translation stage controlled by a stepper motor and a piezoelectric trans- 
ducer. The output power from the mode-locked laser is split and fed to three erbium-doped 
fiber amplifiers (EDFAs). These are used to generate three phase-coherent optical frequency 
combs whose spectral properties can be independently optimized to perform different functions. 
The output from the first EDFA is broadened using a nonlinear fibre to span the wavelength 
range from approximately 1000 nm to 2100 nm. This provides the octave-spanning spectrum 
required for detection of the carrier-envelope offset frequency /o using the self-referencing tech- 
nique [17|18j . With proper adjustment of the polarisation controllers a beat signal with a signal 
to noise ratio SNR=: 40 dB in a resolution bandwidth of 100 kHz is achieved. The electronics 
for the stabilization of the carrier offset frequency comprises a photo detector, a tracking os- 
cillator, and a digital phase-locked loop (PLL) with adjustable gain and bandwidth. The offset 
frequency is stabilized by feedback to the pump laser diode current. The critical frequency to 
be measured is the pulse repetition frequency frep since the optical frequency is measured as 
a very high multiple of this pulse repetition frequency. To overcome noise limitations due to 
the locking electronics and to enhance the resolution of the counting system, we detect frep at 
a high harmonic using a fast InGaAs photodiode with a bandwidth of 25 GHz. Locking of the 
repetition frequency is provided by an analogue phase-locked loop comparing the 90**^ harmonic 
of frep with a microwave reference and controlling the cavity length. The subsequent use of a 
harmonic tracking filter further enhances the short term resolution of our counting system. For 
this purpose, the 9 GHz beat signal is down-converted with a low noise 9 GHz signal synthesized 
from the microwave frequency reference (CSO / hydrogen maser referenced to a Cs-fountain 
clock). The difference frequency is again multiplied by 128, reducing frequency counter digiti- 
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zation errors to below the level of the noise of the microwave reference. Thereby, the frequency 
at which frep is effectively measured is 1.15 THz. A second EDFA generates high power radi- 
ation which is frequency-doubled using a PPLN crystal, generating a narrow-band frequency 
comb around 780 nm. This is subsequently broadened in a nonlinear fiber to generate a comb 
spanning the range 600-900 nm. Light of the spectroscopy laser at 698 nm is superimposed on 
the output of the frequency comb using a beam splitter. To assure proper mode matching of 
the beams and proper polarization adjustment one output of the beam splitter is launched into 
a single mode fiber and detected with a DC photodetector. The other output is dispersed by 
means of a diffraction grating. A subsequent pinhole placed in front of the photodetector then 
selects a narrow range of the spectrum at 698 nm and improves the signal to noise ratio of the 
observed heterodyne beat. For the heterodyne beat signal with the Sr clock laser a SNR of BO- 
SS dB in a bandwidth of 100 kHz was achieved. Again, a tracking oscillator is used for optimal 
filtering and conditioning of the heterodyne signal. All beat frequencies and relevant AOM- 
frequencies were counted using totalizing counters with no dead-time; the counters correspond 
to 7T-estimators for the calculation of the standard Allan variance. 



3 First order Zeeman effect 

In presence of a magnetic field, both clock levels, which have a total momentum F = 9/2, are 
split into 10 Zeeman sub-states. The linear shift Az of a sub-state due to a magnetic field B is 

Az = ynpgElJ-BB/h, (1) 

where mp is the Zeeman sub-state (here 9/2 or —9/2), gp the Lande factor of the considered 
state, fjLB the Bohr magneton, and h the Planck constant. Using the differential g- factor between 
■^Pq and ^Sq reported in Ref. [20]: Ag = 7.77(3) x 10~^, we can determine the magnetic field 
by measuring two symmetrical resonances. Fig. [3] shows the typical resonances observed with 
a magnetic field B = 87^T for both rap — ±9/2 sub-states. The linewidth is of the order 
of 30 Hz, essentially limited by Fourier broadening hence facilitating the lock of the frequency 
to each resonance. This linewidth corresponds to an atomic quality factor Q = 1.4 x 10^^. At 
this magnetic field, two successive 7r-transitions are separated by 96 Hz, which is high enough 
to entirely resolve the Zeeman sub-structure with that type of resonance and to limit possible 
line-pulling effects to below 10~^^ (see section H3| . 

The magnetic field used for pumping and detecting the atoms is provided by two coils in 
Helmoltz configuration to produce a homogeneous field at the center of the trap. They are 
fed by a fast computer-controlled power supply to reach the desired value in a few ms. This 
setup requires to accurately characterize the stability of the magnetic field, as the residual 
magnetic field fluctuations are a possible issue for the clock accuracy and stability. The Zeeman 
effect can provide a precise measurement of this field and its calibration when we measure 
the clock transition for two symmetrical transitions. When probing the two transitions for the 
mp = ±9/2 sub-states, the difference Av between the two frequencies can be related to the 
magnetic field using Eq.[T] Av = ^Ag^sB /h. 

To evaluate the stability of the magnetic field, we chose a particular set of parameters 
[B = 87 /^T and a modulation depth of the numerical servo-loop of 10 Hz) and repeated a large 
number of times the corresponding time sequence as described in section l^T^ The measured 
magnetic field is averaged over 64 cycles. We then concatenated all the averaged data and 
calculated the Allan standard deviation to determine the long term stability of the magnetic 
field. The result is plotted on Fig.[31 The deviation is 10"^ in fractional units at 32 s, and going 
down following a r^^/^ law for longer times. The deviation for long times is below 10^"^. This 
measurement is totally dominated by the frequency noise of the Sr clock and no fluctuations of 
the field itself are visible at the present level of resolution. For a magnetic field of 87 /^T, this 
represents a control of the magnetic field at the sub-/LtT level over long timescales. 
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Fig. 3. Experimental resonances observed for the mp ~ 9/2 mp ~ 9/2 (left) and mp ~ —9/2 ^ 
mp = —9/2 (right) transitions for a magnetic field B = 87 /iT and an interrogation time of 20 ms. 
The lines are gaussian fits to the data. The asymmetry between both resonances results from the 
imperfection of the optical pumping. 
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Fig. 4. Allan standard deviation of the magnetic field. The line is a t~^^^ fit to the data. 
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Fig. 5. (a) Clock frequency as a function of the applied magnetic field. The line represents a fit of the 
experimental data by a quadratic law with one adjustable parameter : the frequency offset. The linear 
term was set to and the quadratic term to its theoretical value, (b) Clock frequency after correction 
for the second order Zeeman efi'ect. The line is the average of the data. 



4 Frequency accuracy 

4.1 Second order Zeeman effect 

The clock has been operated with different values of the bias magnetic field up to 0.6 mT. As 
explained before, our method of interrogation makes the measurements independent on the first 
order Zeeman effect. On the other hand, both resonances are shifted by the same quadratic 
Zeeman shift which has to be evaluated. From the calibration of the magnetic field, we can 
evaluate the dependence of the transition frequency as a function of this field. The results are 
plotted on Fig.O The line plotted on the graph represents the expected quadratic dependence 
of —23.3 Hz/mT^ [21] where we adjusted only the frequency offset to fit the data. The statistical 
uncertainty on this fit is 0.2 Hz, and there is no indication for a residual first order effect to 
within less than IHz at 0.6 mT. At B = 87 /xT, where most of the measurements were done, 
the correction due to the quadratic Zeeman effect is 0.1 Hz only. Conversely, an experimental 
value for the quadratic Zeeman effect coefficient can be derived from the data plotted in Fig.E] 
with a 7% uncertainty. We find -24.9(1.7) Hz/mT^, which is in agreement with theory. 



4.2 Residual lattice light shift 

The clock frequency as a function of the trapping depth is plotted on Fig. [SI Measurements have 
been done with depths ranging from 50 to 500 corresponding to an individual light shift of 
both clock levels levels up to 1.8 MHz. Over this range, the scatter of points is less than 2 Hz 
and the statistical uncertainty of each point lower than 1 Hz. The control of this effect has been 
evaluated by fitting the data with a line. The slope represents a shift of 0.5(5) Hz at 500£'r-- 
The differential shift between both clock states is therefore controlled at a level of 3 x 10~^. In 
ultimate operating conditions of the clock, a trapping depth of 10 -E^ is theoretically sufficient to 
cancel the motional effects down to below 10~^^ [TJ. The light shift corresponding to this depth 
is 36 kHz for both level, or 8 x 10~^^ in fractional units. The kind of control demonstrated here 
would correspond, in these ultimate conditions, to a residual light shift below 2 x 10^^^. 



4.3 Uncertainty budget 

Other systematic effects have been evaluated and included in the accuracy budget listed in 
TablelH 
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Fig. 6. Clock frequency as a function of the dipole trap depth in terms of recoil energies. On the upper 
scale is the corresponding light shift of the clock levels. The line is a linear fit to the data. The value 
of the light shift due to the trap at 500 Er is only 0.5(0.5) Hz. 



The line pulling by neighbouring transitions has been carefully evaluated. Two types of 
transitions should be considered: transverse motional sidebands and transitions between the 
various Zeeman states of the atoms. 

Transverse motional sidebands can be excited by the transverse fc-content of the probe 
laser. For a lattice depth of 100 S^, the transverse oscillation frequency is about 150 Hz. Both 
diffraction and misalignement are below 1 mrad here so that the transverse dynamics is deeply 
in the Lamb-Dicke regime and the height of transverse sidebands are expected to be at most 
5 X 10""^ of the carrier (experimentally, they do not emerge from the noise of the measurements). 
The corresponding line pulling is therefore below 0.4 Hz. This is confirmed by the absence of 
pathological dependence of the clock frequency as a function of the lattice depth (Fig. [6]). 

Unwanted Zeeman transitions result from the imperfection of the optical pumping process 
and of the polarization of the probe laser. In standard configuration the only visible stray 
resonance is the mp = ±7/2 — mp = ±7/2 transition, with a height that is about half of the 
one of the mp = ±9/2 — mF = ±9/2 resonance. It is difficult to set a realistic theoretical upper 
limit on this effect since we have no direct access to the level of coherence between the various 
mp states, nor on the degree of polarization of the probe laser. Experimentally however, several 
parameters can be varied to test the effect. The magnetic field dependence shown in Fig.Oshows 
no deviation from the expected law to within the error bars. On the other hand, measurements 
performed with various depths of the servo-loop modulation also show no differences to within 
0.5 Hz. Finally, we operated the clock with a probe laser polarization orthogonal to the bias field 
and using the rap — ±9/2 — rap = ±7/2 transitions as clock resonances. The clock frequency 
in this configuration is found 0.2(5) Hz away from the frequency in the standard configuration. 
These measurements also test for a possible line pulling from higher order stray resonances 
involving both a change of the transverse motion and of the internal Zeeman state which can 
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Fig. 7. Two series of measurements performed with different clock parameters (see text). The series 
plotted on the right hand side of the figure clearly exhibits a scatter of points that is incompatible 
with the individual statistical error bars of the measurements. Its reduced 'x is 4.3. The other series 
behaves normally and is shown for reference. 



fall very close to the main resonances. Altogether we estimate a 0.5 Hz uncertainty on these 
line pulling effects. 

Incidently, the combination of the measurements using a and tt transitions allows the deriva- 
tion of the differential Landc factor between both clock states [20i. We find g(^PQ) — gCSo) — 
7.90(7) X 10~^, a value that differs from the one reported in Ref. 20] by twice the combined 
1-sigma uncertainty. 

The light shift due to the probe laser is essentially due to the off-resonant coupling of the 
^Pq state with the '^S'l state. The typical light intensity used for the clock evaluation was of a 
few mW/cm^. By varying this intensity by a factor up to 2, no visible effect has been observed 
to within the uncertainty. A more precise evaluation was carried out using the bosonic isotope 
®*Sr [22]. The measured light shift for an intensity of 6W/cm^ was in this case of —78(11) Hz. 
The corresponding effect for our current setup, where the probe power is 3 orders of magnitude 
smaller, is about 0.1 Hz with an uncertainty in the 10~^Hz range. 

The Blackbody radiation shift is derived from temperature measurements of the vacuum 
chamber using two Pt resistors placed on opposite sides of the apparatus and using the accurate 
theoretical calculation reported in Ref. [23 . The Blackbody radiation shift in our operating 
conditions is 2.39(10) Hz. 

Finally a 1 Hz uncertainty is attributed to an effect that has not been clearly identified. After 
having varied all the parameters necessary for estimating the systematic effects, we decided to 
check the overall consistency by performing three series of measurements with fixed parameters. 
The two first ones were performed with a bias field of 87 /iT and servo loop modulation depths of 
7 and 10 Hz respectively. The third one with a larger field of 140 /iT and a modulation depth of 
7 Hz. The results of series 2 and 3 are shown in Fig. [71 where the error bars include the statistical 
uncertainty of each measurement only. The scatter of points of series 3 is clearly incompatible 
with the individual error bars (the reduced of this distribution is 4.3). In addition, its average 
value is 1.5 Hz away. Having not clearly identified the reason for this behaviour (one possibility 
could be a problem in the injection locking of one of the slave lasers at 698 nm), we decided 
to keep this series of data. We also cannot ensure that the effect is not present (though at a 
smaller level) in the other measurements and decided to attribute an uncertainty of 1 Hz to this 
effect. 

Taking into account these systematic effects, the averaged clock frequency is determined 
to be i^ciock =429 228 004 229 873.6(1.1) Hz. The global uncertainty, 2.6 x 10"^^ in fractional 
units, corresponds to the quadratic sum of all the uncertainties of the systematic effects listed 
in Table[T] The statistical uncertainty is at the level of 0.1 Hz. 
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Table 1. Uncertainty budget. 



Effect 


Correction (Hz) 


Uncertainty (Hz) 


Fractiona,! uncertainty (xlO 


Zeeman 


0.1 


0.1 


0.2 


Probe laser Stark shift 


0.1 


< 0.1 


< 0.1 


Lattice AC Stark shift (100 Br) 





0.2 


0.4 


Lattice 2nd order Stark shift (100 Br) 





0.1 


0.2 


Line pulling (transverse sidebands) 





0.5 


1.1 


Cold collisions 





0.1 


0.2 


Blackbody radiation shift 


2.39 


0.1 


0.1 


See text 





1 


2.3 


Fountain accuracy 





0.2 


0.4 


Total 


2.59 


1.1 


2.6 



5 Conclusion 

We have reported here a new measurement of the frequency of the ^ '^Pq transition of 
®^Sr with an uncertainty of 1.1 Hz or 2.6 x 10""'^^ in fractional units. The result is in excellent 
agreement with the values reported by the JILA group with a similar uncertainty [I] and by the 
Tokyo group with a 4 Hz error bar [5] . Obtained in independent experiments with significant 
differences in their implementation, this multiple redundancy strengthens the obtained results 
and further confirms the possibility to build high accuracy clocks with cold atoms confined in 
an optical lattice. It also further assesses this transition as a possible candidate for a future 
redefinition of the second. 

SYRTE is Unite Associee au CNRS (UMR 8630) and a member of IFRAF. This work is 
supported by CNES and DGA. PTE acknowledges financial support from the German Science 
foundation through SFB 407. 
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